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Abstract 

■We report experimental upper limits on WIMP-nucleon elastic scattering cross sections from the second science run of 
|ZEPLIN-III at the Boulby Underground Laboratory. A raw fiducial exposure of 1,344 kg-days was accrued over 319 days 
of continuous operation between June 2010 and May 2011. A total of eight events was observed in the signal acceptance 
region in the nuclear recoil energy range 7-29 keV, which is compatible with background expectations. This allows the 
exclusion of the scalar cross-section above 4.8xl0~^ pb near 50 GcV/c^ WIMP mass with 90% confidence. Combined 
with data from the first run, this result improves to 3.9xl0~^ pb. The corresponding WIMP-neutron spin-dependent 
cross-section limit is 8.0 xlO^'^ pb. The ZEPLIN programme reaches thus its conclusion at Boulby, having deployed and 
exploited successfully three liquid xenon experiments of increasing reach. 



JCeywords: ZEPLIN-III, dark matter, WIMPs, liquid xenon detectors 



(N 
> 

ON 

d 



X 



i. Introduction 

] Direct, indirect and accelerator searches for ncutralino 
dark matter are now probing regions of parameter space 
'favoured by minimal supersymmetric (SUSY) extensions 
to the standard model, in particular those constrained at 
|the unification scale and by measurements of cosmologi- 
cal cold dark matter abundance. SUSY is motivated by 
|the need to stabilise the weak scale, but it is remarkably 
■persuasive that i?-parity conserving flavours of the the- 
ory lead to an excellent WIMP dark matter candidate in 
the form of the lightest SUSY particle. However, at a 
|time when no evidence for SUSY has yet emerged at the 
it is worth noting that direct searches, 
such as the one reported here, aim to detect any WIMP, 
not just neutralinos. Such experiments exploit the possi- 
bility that WIMPs may scatter off ordinary baryonic mat- 
ter. The experimental challenge lies in conducting a rare 
event search (<1 evt/kg/yr) whilst maintaining efRcient 
detection of very low energy signatures (few keV). Liquid 
xenon (LXe) is an excellent target material for intermedi- 
ate mass WIMPs due to its high atomic mass and sensitiv- 
ity in two response channels (scintillation and ionisation). 
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Significantly, these allow discrimination between electron 
recoils resulting from radioactivity backgrounds and the 
nuclear recoils expected from WIMP elastic scattering. 

The ZEPLIN-III experiment operated at the Boulby lab- 
oratory (UK) under a rock overburden of 2,850 m water 
equivalent. This two-phase xenon emission detector mea- 
sures both scintillation and ionisation responses from par- 
ticle interactions in its 12-kg LXe target. Approximately 
half of this mass forms a 'fiducial' region with well un- 
derstood performance and backgrounds. The ionisation 
released at an interaction site is drifted upward and emit- 
ted into a thin (few mm) vapour phase above the liquid, 
where it is converted into an optical signal via electrolu- 
minescence. This is achieved with a strong electric field 
of 3-4 kV/cm in the liquid phase (approximately twice as 
strong in the gas). An array of 31 photomultiplier tubes 
(PMTs) is located within the cold liquid and views the 
36.5-mm thick active region above it. The array responds 
to the prompt scintillation and the delayed electrolumines- 
cence signals (termed SI and S2, respectively). A detailed 
description of the detector design and construction can be 
found in Refs. 0, Q- This time projection chamber config- 
uration allows very good position reconstruction in three 
dimensions, as well as electron/nuclear recoil discrimina- 
tion which is critical for WIMP searches [i,[i,[l^. 
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The first science run (FSR) of tlie experiment in 2008 
placed very competitive upper limits on the WIMP- 
nucleon scattering cross sections in several interaction 
models ll|, [ij, ll3| . The FSR sensitivity was limited by 
background originating from PMT 7-rays. In particular, 
the most challenging event topology in ZEPLIN-III comes 
from multiple-scintillation single-ionisation (MSSI) events, 
whereby a single interaction vertex in the active region 
(producing both SI and S2) is accompanied by one or 
more scatters in a region yielding no charge (SI only). As 
the scintillation responses are effectively time-coincident, 
S2/S1 ratios reconstructed for these 7-ray events are essen- 
tially lower than typical for single-scatter electron recoils, 
and they can leak down to the nuclear recoil acceptance 
region. Nevertheless, we were able to achieve an average 
electron/nuclear recoil discrimination power of 7,800:1 in 
the 2-16 keVee WIMP search region, which is the best 
reported for a LXe detector (hereafter. 'keVee' represents 
the electron-equivalent energy as calibrated by 122 keV 
^^Co 7-rays and 'kcVr' denotes nuclear recoil energy). 

Two upgrades of the experiment had been planned from 
its inception. The first was the replacement of the PMT 
array, which dominated the 7-ray and neutron background 
budgets in the FSR by a large factor. A new PMT model 
(ETEL D766Q Q) was developed in collaboration with 
the manufacturers, which delivered a 40- fold improvement 
in 7-ray activity per unit. This allowed an 18-fold reduc- 
tion in overall electron recoil background at low energies 
relative to the FSR. In Ref. [l^ we analysed the radioactiv- 
ity backgrounds affecting the experiment in the second run 
and showed that they were predicted with good precision. 
Unfortunately, the optical and electrical performances of 
the new PMTs were substantially poorer than those of 
the previous tubes; worse still, the dispersion of gains and 
quantum efhciencies (QE) posed very considerable prob- 
lems to data analysis. Obtaining a set of working PMTs 
and coping with this variability across the array became 
the main challenge of the second run. 

The second major upgrade was the addition of an anti- 
coincidence veto system, which was retrofitted around the 
ZEPLIN-III target [l3,[l3- This 52-module plastic scintil- 
lator detector envelops a Gd-loaded polypropylene shield 
(wStt coverage) which provides moderation and radiative 
capture of internal neutrons; the mean capture time is 
10.7±0.5 /iS. The assembly fits inside the lead shield used 
in the FSR. The veto provides 60% tagging efficiency for 
internal neutrons, mostly derived from a 0.2-70 ^s delayed 
coincidence window from the ZEPLIN-III trigger point 
(58%); elastic recoils in the plastic scintillator within a 
narrow prompt window (±0.2 /is) make up the remaining 
neutron efficiency. The tagging efficiency for 7-rays in the 
prompt window was 28%. Adding to its background rejec- 
tion capability, this tonne-scale detector provides a useful 
source of diagnostic for the radiation environment around 
the instrument. 

Minor upgrades were also implemented to aid with cal- 
ibration of ZEPLIN-III. The radioactive source delivery 



system was fully automated. A copper structure ('phan- 
tom' grid) was installed above the anode mirror; this cast 
a shadow from ^''Co 7-rays onto the LXe surface, thus pro- 
viding calibration of the position reconstruction algorithm. 
A new fibre-coupled LED light gun helped with calibra- 
tion of the PMT single photoelectron responses. All of the 
above upgrades were manufactured from low background 
components. 

2. Second science run 

In the second science run (SSR), WIMP-search data 
were acquired over 319 days between 24*'' Jun 2010 and 
7*^* May 2011, giving a fiducial exposure of 1,344 kg-days. 
A 20% reduction from the FSR fiducial mass to 5.1 kg was 
motivated by the poor performance of peripheral PMTs. 
Even so, the SSR accumulated 3 times more exposure, 
and achieved the longest continuous WIMP run of a xenon 
detector to date. A daily operational duty cycle of 96% 
was achieved consistently, with 1 hr per day reserved for 
^^Co calibration and cryogen re-filling; these tasks were 
automated and controlled remotely. On a weekly basis, 
the system levelling was adjusted (to compensate for lo- 
cal geological movement) and the target and veto PMTs 
were calibrated with their respective fibre-coupled LED 
systems. The gas phase was ^3.5 mm thick and kept at 
1.6 bar, with 13 mbar rms variability within the dataset. 
The electric field in the liquid was 3.4 kV/cm. A free 
electron lifetime of 14 ^s was achieved by prior purifica- 
tion with a hot getter (comparable to the 13-/xs drift time 
for the deepest fiducial interactions). As observed in the 
FSR, the lifetime increased steadily without external pu- 
rification during the run, and eventually reached ^45 /is. 

Data acquisition procedures were described with the 
FSR results [ll[ and we provide only a brief summary 
here. The ZEPLIN-III trigger is derived from the 31-PMT 
sum signal; the SSR hardware threshold corresponded to 
the electroluminescence of si5 electrons; this translates to 
«10 ionisation electrons in the liquid (average over fidu- 
cial interactions in the dataset) when the finite electron 
lifetime and 66% emission probability at the surface are 
accounted for. Single emitted electrons generate a mean 
response of 11.8±0.4 photoelectrons (phe); a study of the 
single electron signature can be found in Ref. [l^l ■ This S2- 
derived signal also triggers the veto data acquisition. The 
two data streams are recorded and reduced separately, and 
synchronised offline. The ZEPLIN-III waveforms are digi- 
tised at 500 MS/s for ±18 /is around the trigger point. Key 
detector and environmental data from the 'slow control' 
acquisition system were embedded into the main data. 

2.1. Data analysis 

Pulse-finding and parametrisation of the waveforms 
were carried out by ZE3RA [l^ , our event reduction soft- 
ware. The reduced data were then searched for single scat- 
ter events (single SI and S2 pulses) which were retained for 
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further analysis. Several corrections were applied, mainly 
to the S2 response, based either on the slow-control infor- 
mation embedded with each event, or on historical trends 
derived from the daily calibration. The mean electron 
lifetime correction was 37%, averaged over the fiducial 
dataset; also corrected were the electronics gain drift (4.7% 
rms), detector tilt (1.9%) and pressure variations (1.1%). 

A vertex reconstruction algorithm estimated the energy 
(E) and position {x,y) of the interactions, with a maxi- 
mum likelihood fit to the SI response and a least squares 
fit to the S2 channel [2^. The algorithm fits to 31 empir- 
ical response functions simultaneously. These model the 
response of a PMT as a function of distance to its axis, 
and are derived from calibration data in an iterative pro- 
cedure which also 'flat-fields' the array. The spatial res- 
olution thus achieved was 13 mm in SI and 1.6 mm in 
S2 (FWHM) for the inner 100 mm radius. The iE,x,y) 
parameter space was navigated with a simplex method to 
obtain the best fit point for each event; spatial maps of 
likelihood and could also be produced to help identify 
multiple vertices. 

Finally, fiducialisation and quality cuts were applied to 
the data. A central region of 140 mm radius contain- 
ing 5.1 kg of LXe was retained. The quality cuts elimi- 
nate likely MSSI events and outliers in several parameters 
(e.g. goodness-of-fit of SI and S2 reconstructions). An SI 
pulse timing cut was also used to exploit modest pulse 
shape discrimination at the higher energies (2l| . The in- 
efficiencies incurred from these cuts are described in the 
next section. Unsurprisingly, they are more severe than 
had been required for FSR data. 

A blind analysis was conducted during the initial opti- 
misation of the quality cuts and definition of the WIMP ac- 
ceptance region: low energy events in and around the sig- 
nal region were kept hidden from visualisation and reduc- 
tion using a 'blindness manager' implemented in ZE3RA. 
Vetoed events (mostly 7-rays) were excluded from this 
list, thus providing a (practically) signal-free population 
of background events. This blind analysis was not pur- 
sued to its ultimate conclusion but the acceptance region 
(which had been defined blindly) was retained in order to 
avoid major bias in estimating a potential signal. 

2.2. Calibration and efficiencies 

Daily calibration with ^''Co 7-rays (mostly 122 kcV) de- 
fined the SI and S2 energy scales and monitored their 
stability throughout the run; it also provided a regular 
measurement of the detector tilt (through the spatial de- 
pendence of the S2 pulse width) and of the free electron 
lifetime (by S2/S1 ratio with interaction depth). The lat- 
ter parameter was confirmed independently with a new 
method based on single electron counting operating on the 
actual science data waveforms [l8l |. 

In the SSR the light yield of the chamber (fiducial aver- 
age) was 1.3 phc/keVec, down from 1.8 phc/keVee in the 
FSR. The energy resolution was similarly affected; a linear 
combination of SI and S2 responses which exploits their 
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Figure 1: Discrimination parameter as a function of Sl-dcrived en- 
ergy for Am-Bc neutron calibration. The recoil population median 
and it2-(T lines are shown (in blue), along with the 2-12 kcVce search 
box (dashed green). 



microscopic anti-correlation yielded 16.4% (FWHM) for 
the new array, against a very impressive 8.4% achieved for 
FSR data using the same analysis codes. In thin LXe tar- 
gets like ZEPLIN-III individual PMTs collect a large frac- 
tion of the scintillation light of an event; the uniformity of 
response in the array becomes therefore more critical than 
in high-reflectance chambers with deeper geometries. In 
spite of this, a spatial resolution of 1.6 mm (FWHM) in 
the horizontal plane was achieved in the SSR, measured 
from the shadow pattern cast by the new phantom grid 
for "Co 7-rays [2^. 

The electron-recoil population was calibrated with a 
4.6 kBq ^"^^Cs source located above the instrument, pro- 
ducing a rate of 150 c/s in the detector. The number 
of low-energy events thus obtained was similar to that of 
background electron recoils in the search data. 

Three neutron calibrations took place during SSR data- 
taking, totalling 10 hours of exposure to an Am-Be 
(a, n) source emitting 1,321±14 n/s [22]; absolute differ- 
ential rates of nuclear recoils from elastic scatters agreed 
within statistical errors for these datasets. The delayed- 
coincidence veto efficiency for such events was 58% inde- 
pendently of recoil energy [l^l (the delayed detection is 
provided by ^^^Gd 7-rays rather than by the neutrons 
directly). One such run is depicted in Figure [TJ The 
recoil band populated by neutron elastic scattering was 
parametrised by Gaussian fitting to the discrimination 
parameter logjQ(S2/Sl) in 1 keVee (SI) bins. Energy- 
dependent mean (/i) and width (cr) parameters were used 
to define a signal acceptance region in the range 2-12 keVee 
and including approximately the lower half of signal accep- 
tance as indicated in the same figure: from 2.3% (/i — 2(t) 
to 45% (/X — 0.126cr). This region was defined after a grad- 
ual unblinding of surrounding data in the WIMP search 
dataset, with outer regions progressively discarded and 
then opened for analysis as 'sidebands'. 
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Figure 2: Nuclear recoil efficiencies. Above: relative scintillation 
efficiency for nuclear recoils in LXe at zero electric field; the thick 
blue curve (adopted) combines the measurements and uncertainties 
reported in Ref. [23[ for FSR data (68% confidence region, in red) and 
SSR data (green). Below; recoil detection efficiencies; rj subsumes a 
flat component »yo=39.8% (dominated by signal acceptance in S2/S1) 
and three energy-dependent factors: rji is the SI detection efficiency 
from the 3-fold PMT coincidence required in software; J72 is incurred 
from quality cuts; rja relates to the timing cut on SI pulses. 



New measurements of the relative scintillation efficiency 
for nuclear recoils in LXe (Lpjf) were derived from FSR 
and SSR calibration data [231 . The decreasing Lg// curves 
obtained therein were mutually consistent and agree, 
within experimental errors, with recent neutron beam 
25| . For the purpose of converting between 



data [24 

electron and nuclear recoil energies we combined the two 
ZEPLlN-111 measurements into the curve shown in Fig- 
ure O (upper panel) . The SSR WIMP search region corre- 
sponds to 7.4-29 keVr. 

The nuclear recoil detection efficiency comprises sev- 
eral components. A set of constant factors combine to 
7;o=39.8%, including: DAQ livetime fraction (99.2%), 
waveform quality cuts (98.3%), cuts on robustness of 
pulse parametrisation (96.9%), veto random coincidences 
(99.6% and 99.0% for the prompt and delayed windows, 
respectively) and the above signal acceptance fraction 
(42.7%). Three additional curves, shown in Figure[5] (lower 
panel), describe the detection efficiency for SI pulses (3- 
fold requirement in software), the efficiency of the data 
quality cuts and the additional timing cut on SI signals 
for pulse shape discrimination. The effective exposure for 
a 50 GeV/c2 WIMP was 251 kg-days. 

2.3. Experiment backgrounds 

Nuclear recoils from neutron single elastic scattering 
pose an irreducible background to direct WIMP searches. 
In ZEPLIN-III, detected radioactivity neutrons are most 
likely to arise in the ceramic feedthroughs, the PMTs 



and the laboratory rock; muon-induced neutrons con- 
tribute negligibly. The overall predicted rate is 3.05±0.5 
events/year in 5~50 keVr assuming unity detection ef- 
ficiency, as detailed in Ref. Translated to the 
SSR effective exposure and search region this gives only 
0.06±0.01 events in anti-coincidence with the veto. 

A low-energy electron recoil background rate of 
0.75 evts/kg/day/keVee ('dru') was measured in the fidu- 
cial volume (with no quality cuts applied), which repre- 
sents a 20-fold reduction brought about by the new photo- 
tubes. Monte Carlo predictions from a comprehensive in- 
ventory of background sources, based on component-level 



radio-assays, indicated 0.86±0.05 dru [15|. Background 
electron recoils are promptly tagged in the external veto 
with 28% efficiency below '-^100 keVee. This confirms that 
internal /3~ emitters are insignificant, since the veto effi- 
ciency for coincident 7-rays is practically identical to this 
value. Independently, we measured the *^Kr decay rate at 
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Assuming the discrimination factor achieved in the first 
run, electron-recoil event leakage into the WIMP accep- 
tance region should represent <1 event, in line with the 
aspiration of a background-free second run which moti- 
vated the experiment upgrades. However, analysis of the 
prompt-vetoed 7-rays revealed that this was unlikely to be 
the case. Background predictions (discussed below) based 
on extrapolation of the electron recoil population into the 
signal box and from low-energy ^'^''Cs data indicated that 
a handful of unvetoed events (7-9) were expected, con- 
firming a loss of discrimination power. 

3. WIMP-search results 

Upon unblinding, 12 events were observed in the ac- 
ceptance region. A detailed waveform inspection revealed 
larger than expected cross-talk artifacts due to the poor 
electrical performance of the SSR PMTs: gains differed 
by as much as 100 times within the array, and the extra 
amplification required in some channels exposed contami- 
nation from the higher gain channels (the array is powered 
by a single HV supply using common dynode electrodes 
internally). Tighter cuts on the goodness-of-fit of the re- 
constructed vertex were required to deal with this issue 
(Figure [2] already reflects this) . The result of the final 
(non-blind) analysis is shown in Figure [31 Eight events re- 
mained in the box. The event reconstructed at 3.2 keVee 
in SI and 1.1 keVee in S2 appears far below the mean 
log]^Q(S2/Sl) for typical electron recoils with that SI sig- 
nal, but this is not necessarily anomalous: a median S2 
signal of 1.1 kcVee corresponds to electron recoils with 
only 0.6 kcVee in SI (cf. yellow line in Figure [3]). Given 
1.3 phc/keVee scintillation yield, 0.8% of all events with 
that SI expected would generate >4 phe due to Poisson 
fluctuation, thereby producing the observed log]^Q(S2/Sl) 
ratio or lower. Considering the detected rate of electron 
recoils, one such low-lying event is consistent with back- 
ground. 
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Table 1: Observations (riobs), background estimates {(11,1,2) and limits on the signal expectation {fJ.s) for the first (FSR) and second (SSR) 
runs of ZEPLIN-III. Fiducial and net effective exposures (50 GcY/c^ WIMP) arc presented along with relevant signal acceptance parameters. 
Electron recoil background expectations are estimated from Skew-Gaussian (SG) fits to data above the search region in 2 keVee bins and also 
from Cs calibration. The 90% CL limit on the number of signal events in each run is derived with the Profile Likelihood Ratio method. 
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kcVr 
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SG fit 
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90% CL 


FSR 


437.0 


7-35 


29-50% 
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0.5±0.3 


5.2±3.1 




<4.2 




(107.3) 




2-29% 
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0.7±0.3 


1.5±1.7 






SSR 


1,343.8 


7-29 


24-45% 


7 


0.03±0.005 


5.5±2.2 


8.3±2.9 


<5.1 




(251.0) 




2-24% 


1 


0.03±0.005 


1.0±1.2 
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No delayed coincidences were recorded below the nuclear 
recoil median and those registered above it are statisti- 
cally consistent with random coincidences. This allows 
us to set an upper limit of 0.75 neutron events for the 
search region (90% CL), confirming the successful mitiga- 
tion of the neutron background. Two methods were used 
to predict counts from electron recoil backgrounds; both 
are in good agreement with the observation, suggesting 
no significant signal. Binned Skew-Gaussian (SG) fits to 
the log]^g(S2/Sl) parameter above the search region, as in 
the FSR analysis [ll|, predict a total of 6.5±3.4 events 
(cf. 8 observed). ^'^^Cs calibration data provided indepen- 
dent confirmation of this background within the available 
accuracy, with no contamination from signal and no as- 
sumption of a functional dependence for electron recoils; 
this yielded 9 events for 96% equivalent exposure. 

To account for non-uniform background in logjQ(S2/Sl) 
we partitioned the search region in two; this was done ex- 
plicitly to maximise sensitivity given the SG background 
distributions. Note that the contents of the search re- 
gion were not considered in this optimisation. To achieve 
this the partition was incremented systematically and the 
background distribution was sampled by Monte Carlo un- 
der the null hypothesis; the best sensitivity, calculated by 
the Feldman-Cousins (FC) method [l^, is reached at 24% 
acceptance. The observations in the resulting bins were 
seven and one events, as shown in Figure [3] (lower). Ta- 
ble 1 confronts observations and background estimates for 
the partitioned search region. 

Reanalysis of FSR data with updated SSR software, mo- 
tivated by the new L^ff measurements of Ref. [2^ , revealed 
five events in the FSR search region (which we kept un- 
changed from the original analysis). Background expec- 
tations from SG fitting to these data are also presented 
in Table 1. In this instance the ^'^^Cs calibration over- 
predicted the observations very significantly, which may be 
caused by rate-dependent photocathode charging effects in 
this (FSR) run [27[ or other systematic differences between 
the two calibrations. We consider, however, that the SG 
function provides an appropriate description of our back- 
ground as validated by the SSR ^"^^Cs calibration. Once 
the non-negligible neutron contribution is included, the 
FSR prediction is 7.9±4.8 (cf. 5 observed). The same 
data-blind optimisation procedure led to a partition at 



29%; in this run one event was just below that division 
in logjQ(S2/Sl), at 28.8% acceptance. 

3.1. Signal inference 

A confidence interval for the signal expectation (/is) in 
the SSR acceptance region was obtained with a Profile 
Likelihood Ratio (PLR) method which accounts for the 
uncertainty in the background predictions (see, e.g. [28 . 
[2^). We included estimators of the nuisance parameter 
fib2 from both the SG fits and the ^'^^Cs calibration; the 
latter is Poisson distributed and the SG predictions are 
treated as Gaussian truncated at zero. The distribution of 
profile likelihood ratio was determined by Monte Carlo to 
ensure correct statistical coverage. The double-sided 90% 
CL interval for fig was — 5.1 events. The same procedure 
applied to FSR data (without a ^'^''Cs prediction in this 
instance) yielded — 4.2 events with 90% confidence. 

These results are quite robust with respect to the partic- 
ular choice of binning and test statistic. For example, we 
conducted two-bin FC calculations with background pre- 
dictions capped at observation when fih > Uohg (so as not 
to benefit from downward fluctuations of background) and 
using a 10% acceptance upper bin (as in the original FSR 
analysis [HI). This yielded the same 4.2 events at 90% 
CL for the first run and 4.8 for the SSR. These are reas- 
suringly close to the PLR results in Table 1, which treat 
background uncertainties more formally. 

Experimental upper limits on the scalar WIMP-nucleon 
clastic cross section are shown in Figure HI calculated with 
the standard galactic halo model {po~0.3 GeV/c^/cm'^, 
110=220 km/s, Wesc=544 km/s, and km/s) and the 

Helm form factor [s^l parametrised in Ref. [3l|. A mini- 
mum cross-section limit of 8.4xl0~® pb (90% CL, double 
sided) is reached at 55 GeV/c^ WIMP mass for the FSR 
(similar to the original res ult jlll| and slightly lower than 
the result reported in Ref. l23[ from a single-sided 'maxi- 
mum patch' test statistic [32|). The minimum of the SSR 
curve is 4.8xl0~® pb, reached at 51 GeV/c^. Adopting 
Leff curves at the ztlcr levels from Figure[2]docs not affect 
the result at curve minimum significantly, but the value at 
10 GeV/c^ mass varies in the range (2.2-9.3) x 10~^ pb. 

The combined result for the ZEPLIN-III experiment 
with mean L^ff, obtained from a four-bin PLR calcula- 
tion which returns /is <6.0 events for the aggregate expo- 
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Figure 4: 90%-CL limits on WIMP-nuclcon spin-independent cross 
sections from ZEPLIN-III (FSR, SSR and combined) and from 
XENONIOO [3, XENONIO (low energy analysis [H), CDMS-II [3 
and EDELWEISS-II [H- Previous results from the ZEPLIN pro- 
gramme are also indicated l38l. l39ll . In blue we represent the 3- and 
5-(T DAMA/LIBRA contours (2008 data, no ion channelling ^) in- 
terpreted in Ref. [4lll . The magenta contour is the fit to CoGeNT 
data under the light WIMP hypothesis l42l : in cyan is the 2-a region 
from CRESST-II [H. The crosses are the original SUSY bench- 
mark points [44l . Favoured regions of parameter space from a 2008 
Bayesian analysis in mSUGRA [45ll and the likelihood analysis of 
LHC data within cMSSM 0| are also shown. 



Figure 3: Above; fiducial events in full SSR exposure. Green mark- 
ers label prompt veto coincidences (mostly 7-rays); events in the 
delayed window (mostly accidental coincidences) are in red. There 
are 8 unvetoed events in the WIMP acceptance region. The dashed 
lines show the nuclear and electron recoil band medians (in blue 
and yellow, respectively). Below: Distribution in signal acceptance 
(recoils from WIMPs and neutrons are distributed uniformly in the 
y-axis). There were seven (one) events in the upper (lower) region. 



sure, is also shown in the figure; the curve minimum is 
3.9x10-* pb at 52 GeV/c^. 

Excellent sensitivity to spin-dependent WIMP-neutron 
interactions is afi'orded by the odd-neutron isotopes ^^^Xe 
and ^•^^Xe. The spin-dependent result is calculated as de- 



scribed in Ref. jl2| , accounting for the composition of our 
xenon (depleted in ^'^^Xe) and using Bonn-CD nucleon- 
nucleon potentials. The FSR+SSR combined curve, shown 
in Figure [5j has a minimum of 8.0x10-"^ pb at 50 GeV/c^ 
mass. At the time of writing there is no corresponding 
result from XENONIOO; this is expected to be a few times 
lower. The original XENONIO result [s^ is not shown in 
the figure since this had assumed a constant scintillation 
efficiency and the more favourable (^^2) Bonn A poten- 
tial; a fair comparison would raise this result to above the 
ZEPLIN-III curve. 



4. Conclusion 

In this article we presented experimental upper limits 
on WIMP-nuclcon clastic scattering cross sections from 
the second run of ZEPLIN-III at Boulby. These were de- 
rived from analysis of 1,344 kg -days of fiducial exposure 
acquired between June 2010 and May 2011. A 90% CL 
upper limit on the signal content from this exposure al- 
lows the exclusion of a scalar WIMP scattering cross sec- 
tion above 4.8x10^* pb/nucleon near 50 GeV/c^ mass. 
The combined result for the two runs constrains the scalar 
cross section to 3.9x10^* pb and the WIMP-neutron 
spin-dependent cross section to 8.0x10^'^ pb. Along 
with XENONlOO^p, XENONIO [H, CDMS-II ^ and 
EDELWEISS-II [37[ these results disfavour an interpreta- 
tion of DAMA in terms of nuclear recoils from WIMPs as 
well as recent results from CRESST-II under the canoni- 
cal dark halo Q. Models favoured by CoGeNT [13] are 
harder to rule out completely from our data. 

This second run followed the upgrade of the experiment 
with a new array of purpose-developed, low-radioactivity 
PMTs, a veto detector based on Gd-loaded polypropylene 
and plastic scintillator, and new calibration hardware. The 
automation of the system enabled a 319-day run with ex- 
cellent stability and reduced manpower underground. 
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Figure 5: Limits on spin-dependent WIMP-neutron cross section 
(assuming no proton interaction) from ZEPLIN-III (FSR+SSR. with 
Bonn CD potentials O) as well as CDMS-II (2004-09 data [^,[33 
plus low-energy (LE) analysis gl), ZEPLIN-I [H and ZEPLIN- 
II [soli . Also shown is the 3-o" DAMA evidence region (2008 data, no 
ion channelling |40|| ) interpreted in Ref . |41|| . The green hatched area 
is the tip of the 95% probability region for cMSSM neutralinos [5l|. 



ing designed and constructed around the world with tonne- 
scale fiducial masses, the ZEPLIN programme can claim 
to have pioneered some of the techniques that helped two- 
phase xenon become a leading technology in the race to 
discover WIMPs. 
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The new photomultipliers were a critical item to reach 
the design sensitivity of 1x10^* pb-year. Although their 
radiological background was excellent (35-50 mBq/PMT), 
their performance compromised the experimental sensitiv- 
ity by a factor of ~4 (from reduction in fiducial mass, 
of signal acceptance fraction, of cut efficiency and from 
poorer discrimination). Electron- recoil leakage to below 
the nuclear recoil median had been 7,800:1 in the FSR 
compared to only 280:1 in the SSR; no other significant 
changes were made to the internal hardware. 

The addition of a new veto detector proved very valu- 
able, confirming negligible neutron background and pro- 
viding a useful unbiased sample of 7-ray background for 
open analysis. In our experience, high-efficiency veto sys- 
tems will be extremely important in future experiments, 
from the point of view of background reduction, diagnos- 
tics and enhancing their discovery power [l7l |. 

ZEPLIN-III concludes a successful series of three differ- 
ent LXe-based experiments operated at Boulby since the 
mid 1990s, with progressively stronger electric fields ap- 
plied to their active targets. ZEPLIN-I exploited pulse 
shape discrimination at zero field and produced world- 
leading results [H,!!^. Delivering 1 kV/cm to the liquid 
target, ZEPLIN-II was the first two-phase xenon WIMP 
detector in the world [H^. ZEPLIN-III operated at 
nearly 4 kV/cm in the first run and achieved the best dis- 
crimination of any xenon detector, along with competitive 
WIMP results in general. As several systems are now be- 
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